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^ : Abstract 

^ ■ We report a fully microscopic theory for transconductivity, or, equiva- 



lently, momentum transfer rate, of Coulomb coupled electron systems. We 



o 

O ■ use the Kubo linear response formalism, and our main formal result expresses 



: ^ ■ the transconductivity in terms of two fluctuation diagrams, which are topolog- 

>< 

■ ically related, but not equivalent to, the Aslamazov-Larkin diagrams known 

for superconductivity. Previously reported results are shown to be special 
cases of our general expression; specifically, for constant impurity scattering 
rates, we recover the Boltzmann equation results in the semiclassical clean 
limit, and the memory function results in dirty systems. Furthermore, we 
show that for energy dependent relaxation times, the final result is not ex- 
pressible in terms of standard density-response functions. Other new results 
include: (i) at T = 0, the frequency dependence of the transfer rate is found 
to be proportional to 17 and for frequencies below and above the impurity 
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scattering rate, respectively and (ii) the weak localization correction to the 
transconductivity is given by Sa2i'" oc Sa^^ + 6a22^- 



Typeset using REVT^iX 



I. INTRODUCTION 



Consider two systems containing mobile charge carriers so close to each other that the 
charges in the two respective subsystems feel the Coulomb forces originating from the other 
subsystem, and yet far enough away from each other that direct charge transfer between the 
two subsystems is not possible. Experimental realizations of such systems are, for example. 
Coulomb coupled double quantum well systems,0i arrangements where a 3D system is close 
to a 2D system,! or two nearby quantum wires. A scattering event between a carrier in 
one system and a carrier in the other system leads to momentum transfer between the two 
subsystems. Thus, if a current is driven through one of the systems (henceforth the driven 
system is denoted as layer 1), then an induced current is dragged in the other subsystem 
(layer 2). Alternatively, if no current is allowed to flow in layer 2, a voltage is induced. Due 
to momentum conservation the two particle number currents flow in the same direction. 
Since the mechanism for the Coulomb drag is carrier-carrier scattering the drag current 
is proportional to the square of the effective interaction between the subsystems. The 
available phase space for electron-electron scattering tends to zero at low temperatures, and 
consequently one expects Coulomb drag to decrease with decreasing temperature. At low 
temperatures, the two Pauli factors entering the carrier-carrier scattering rate lead to a T^- 
dependence, and this behavior is approximately seen in experiments.0 Note, however, that 
there are small, but important deviations from the simple T^-law; these deviations have 
been the topic of much recent interest. 001 

The possibility for Coulomb drag was realized already long ago,ill and the recent ex- 
perimental advanceJlHi have brought about a flurry of theoretical works. A number of 
different theoretical approaches has been proposed. These include (i) calculations based 
on the Boltzmann equation,!'! (ii) the memory function approach of Ref. |^, and (iii) the 
momentum balance equation method.EI 

In this paper, we calculate the Coulomb drag between two systems using a fully micro- 
scopic theory based on a linear response formula. The central object to be evaluated is the 



3 



retarded current- current correlation function; since the two involved currents refer to the 
the two different subsystems, we call the result of this calculation trans conductivity. The 
motivation underlying our work is that all previously proposed approaches lack the rigor 
that can be achieved with a formal linear-response calculation. The present method allows 
us to identify the Feynman diagrams that contribute to the transconductivity. Instead of the 
normal conductivity bubble, we find that one must evaluate a fluctuation diagram, which is 
similar, but not identical, to the Aslamazov-Larkin0 diagrams known from superconductiv- 
ity, or the diagrams encountered in connection with the microscopic theory of van der Waals 
interactions. Thus, all the methods developed within the diagrammatic perturbation 
theory are readily applicable, and one can systematically study the effects of higher or- 
der scattering processes, such as vertex corrections or weak localization, or electron-phonon 
interactions, or the effect of magnetic fields. 

Apart from the general formulation for calculating transconductivities, we obtain the 
following explicit results. In the limit of weak impurity scattering we show that the linear- 
response result reduces to the expression obtained with the Boltzmann equation. We also 
study the corrections to the Boltzmann equation formula in the case of stronger impurity 
scattering, i.e. accounting for vertex corrections. Further, we show that the Boltzmann 
equation result, which involves the susceptibility functions of the individual subsystems, 
must be generalized, if one considers energy-dependent scattering rates. We also consider 
weak localization corrections to the transconductivity. All these effects are calculated at 
finite temperature, but for zero external frequency. At T = finite frequency calculations 
become feasible, and we present a general proof that the dc-drag current vanishes in the 
dc limit at zero temperature for an open system. Recently, Rojo and MahanEl used a 
ground state energy argument to calculate the drag current, and found a nonzero result 
at zero temperature. This seems to contradict our results, however there is an important 
distinction: the calculation of Ref. |T3| applies for closed systems, e.g. coupled mesoscopic 
rings. 

This paper is organized as follows. Section || outlines the derivation of the general ex- 
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pression for transconductivity. Section |T| is devoted to impurity scattering: first, we show 
how the well-known Boltzmann result follows from the general formulation, and next, we es- 
tablish a connection to the memory function formulation of Ref. |^. This section is concludes 
with a discussion of energy-dependent scattering rates and weak localization effects. Section 
presents result for T = at finite frequency. Finally, a number of technical details can 
be found in the appendices. 

II. GENERAL FORMULA FOR THE TRANSCONDUCTIVITY 

The previous works have related Coulomb drag to the transresistivity, p2i] our calcula- 
tion, which is based on a Kubo formula leads to the transconductivity, a2i- These are defined 
as 

P21 = with J2 = 0; (la) 
0-21 = 77-, with E2 = 0, (lb) 

where Ei and Jj are, respectively, the electric field and the current density in layer i. These 
two quantities are related via 

-0-21 -(T21 .^x 
P21 = ^ (2) 

0"ll<722 ~ Cri20"21 0"ll<722 

In (|^) the diagonal a's are the individual subsystem conductivities and we note that the 
transconductivity is always much smaller than intralayer conductivity, because it is caused 
by a screened interaction between spatially separated systems {e.g. the data of Ref. |l] gives 
o"2i/o"ii — 10"^). The transresistivity P21 is a more physically relevant quantity than a2i in 
a drag-rate measurement because p2i is directly related to the rate of momentum transfer 
from particles in layer 1 to layer 2, T2'i, without reference to the scattering rates of the 
individual layers; i.e.. 



mi _ 1 idp2/dt)drs.g 



P21 = — 2 — ; — = z , (3) 

^ie^r2i r2i Pi 
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where m is the effective mass, n is the carrier density, p is the momentum per particle, 
{dp/dt)dra.g is the momentum transfer due to interlayer interactions, and the overhne denotes 
an ensemble average. 

The Kubo formula^ expresses the conductivity tensor in terms of the retarded current- 
current correlation function, 

affix - x'- n) = '■(^ - x'; n) + ^5{x - x')S,,S^pp,ix), (4) 

where (throughout we use ^ = 1) 

U^f^^ix -x'-t- t') = -zQit - t'){[jnx,t),j^{x',t')]). (5) 

Here {ij} indicate the subsystem, {a/3} in the superscripts label the Cartesian coordinates, 
Pi{x) is the particle density in subsystem i, and j{x,t) is the particle current operator. 
We have assumed that the subsystems are translationally invariant. Our task consists of 
calculating the transconductivity ■ 

We employ the imaginary-time formalism to evaluate the retarded current- current cor- 
relation function, starting with the (imaginary-)time-ordered correlation function 

U-^ix -x'-r- r') = -(T.{jT(^, (^', r')}). (6) 
The retarded function then follows as 

n^f '"(a; - x'; Q) = lim H^f (a; - x'; (7) 

where 

n^f (a; - x'; iQ^) = f dr e'^""n°f (a; - x'- r) (8) 

n^f (a; -x'-r) = W e~*^-n^f (x - x'- (9) 

P n 

and where (3 = l/ksT. The calculation proceeds by expanding the transconductivity in 
powers of the interaction between the subsystems. The interaction Hamiltonian is given by 

Hi2 = J dri J dr2 Pi(ri) U^iri - 7-3) P2(r'2)- (10) 



Here Uu is the bare Coulomb interaction between the the systems. We note that other 
interaction processes, which couple the charge carriers in the two subsystems, can be treated 
similarly. An important example is the virtual phonon mediated interaction, which may play 
a role in the low temperature behavior of the momentum transfer rate.B^ 

The r-dependence of the current operators in is determined by the full Hamiltonian 
H = Hi + H2 + H12, where Hi are the subsystem Hamiltonians. In order to develop 
a perturbation expansion, we must isolate the ifi2-dependence. Following the standard 
many-body prescription,lii we transform into the interaction representation and obtain 

n^f (a, _ a;', r - r') = -(T.{^(/3)jr(a., t)j^{x', r')}), (11) 

where the carets indicate that the time-dependence is now governed by the individual sub- 
system Hamiltonians, and the operator S{(3) is 

S{(3) = T.{exp[- J^nHuiri)]}. (12) 

As usual, only connected diagrams need to be included. It is now straightforward to expand 
S{P) in powers of Hu, 

S{f3)^i-T,{jynHi2in)} 

+ IMJo"^^^ dT2Hi2{ri)Hi2{r2)} + . . . (13) 

The zeroth order term leads to a vanishing contribution to the transconductivity because 
the two current operators are decoupled and hence commute. In the following sections we 
discuss the higher order terms. 

A. Linear expansion 

The linear order term in H12 leads to the correlation function 

H^f (a; -x\t- r't^ = dr^ j drirfr2(T.{jT(^, n)}) 

X f/i2(ri-r2)(r.{/52(r2,ri)jf(a^',r')}). (14) 
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Use of the continuity equation, iQp + V • j = 0, allows us to eliminate the number density 
operator, and to express the density-current correlators in terms of the subsystem conduc- 
tivities. After some simplification we find (for a translationally invariant impurity averaged 
system where impurity scattering in the two subsystems is uncorrelated^) 

7,0 

This expression is exact, and it can be used to calculate the first order transconductivity 
for any system, once the subsystem conductivities an are knownJll From ( |15|) we also infer 
that the first-order transconductivity vanishes in the dc- limit. 



B. Quadratic expansion 



To evaluate II21 to second order in H12, we substitute the third term of (|I3D on the 
right-hand side of (|lID , and find that the current- current correlation function is given by 

1 ff^ fl^ 



X 



f/i2(ri - r2)t/i2(r; - r^)A?(a;r, nn, r;r2)A^(a;V, ran, r^rs), (16) 



where we have defined the function 

Af (a;r, x'r\ x"r") = -(T.{jr(a;r)p,(a.V)p,(a;V)}). (17) 

Just as in the previous section, we factorized the time-ordered expectation value involving 
two current and four density operators; this step is justified because the two subsystems are 
decoupled after the formal expansion in H12. Due to the assumed translational invariance A 
depends on only two coordinate differences. We define the Fourier transform A{q,q';uj,uj') 
via (z/ = the volume) 

A{xt,x't',x"t") = — J2 ~~ J2 ^iqiix-x'')+iq2\x'-x'')^-iLo„,{T-T'')-iuj,,{T'-T'') 

xA(qi + q2,q2;iuJm + iuJn,iuJn)- (18) 
The final expression for II21 {Q , ifln)^^'' then takes the form 



A"(Q + q, q, iQ^ + i^n, ^tUn)Af (-Q - g, -q, -iQn - i^n, -iuJn)- (19) 

The diagram corresponding to the second order result, Eg . ([T9|) , is shown in Fig. |l|. We also 
display the first order term discussed above. 

Consider next the iun summations. The function A{z + iQn, z) has branch cuts in the 
complex 2;-plane at Im(if2„ + z)=0, or lm(2;)=0, and is analytic elsewhere (see Appendix 0). 
We can therefore perform the zcij„-sum as a contour integral. When we extract the retarded 
part according to Eq. (0), we obtain the result 



S{Vt) = ^ ^ Ai(ifi„ + iu;„,iu;„)A2(-ifi„ - -iun) 



ILUn 

rco 



= Si{n) + S2in), (20a) 

S,in)=V I [nj^icu + n^iu)] Ai(+, -) A2(-, +) (20b) 

/CO fjfj 
— [nB{uj)A,{+, +)A2(-, -) - n^iu + fi)Ai(-, -)A2(+, +)] , (20c) 
-oo ZTTl 

where n-B{uj) is the Bose function. We have suppressed the Cartesian indices and the mo- 
mentum labels, since they can be gleaned from Eq.(|T9|), and have used the notation 

Ai(±, ±) = A{iQn + ii^n — > + uj±i6; iun uj±i5), (21a) 
A2(±, ±) = A[—iVLn — ioJn —>■ — uj±i6, —iuJn —u!±i6) . (21b) 

The functions Aj(+, +) and Aj(— , — ) vanish identically in the dc limit, which is proven in 
Appendix 0. 

Consider next the dc-response, ^ 0, of a uniform system {Q = 0). The dc-limit of 
Si{fl) is simple to evaluate, because the difference of the two Bose functions combined the 
prefactor ^2^^ just gives dujUB^uj). Thus the dc-transconductivity reduces to (reintroducing 

») 

X A"(q, q;uj + i6,uj — i6)A2{—q, —q, —uj — i6, —u + i6) . (22) 
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The actual evaluation of this expression at various levels of approximation forms the main 
task of this paper. In the case of electron-hole systems the overall sign of (^) must be 
changed. 

C. Higher order terms 

The S-matrix expansion (0) can be used to generate higher order terms. To proceed 
systematically one must apply the techniques of the many-body formalism. As usual, the 
most important processes should be identified, and the corresponding diagrams be summed 
to infinite order. This procedure may then lead to an integral equation for the effective 
interaction, e.g. in the ladder approximation one obtains the Bethe-Salpeter equation. We 
do not pursue this line of argument further in this paper, but note that a particularly useful 
resummation can be obtained, if one includes the "bubble" -diagrams (see Fig. Q), which 
leads to an effective screened interaction, Ui2{q) Ui2{q,uj) = f/i2(g)/ei2(g, cu), where the 
dielectric function is given by 

ei2(g,w) = [1 - Ui{q)xi{q,uj)][l - U2{q)x2{q, uj)] - Ui2{qfxi{q,^)X2{q,^), (23) 

where the x's are the usual polarization functions, and the f/j's are the intrasystem Coulomb 
interactions. We observe that an energy-dependent U (g, u) can be used in the above ex- 
pressions, (|15]) and (^2]), for transconductivity with no additional difficulty. Most previous 
wor on drag problems have used ( ^31) (or simplified versions of it). 

III. IMPURITY SCATTERING 

In the previous section we showed that the transconductivity can be expressed in terms 
of the general three-body correlation function A. We will next consider a specific example 
in order to calculate this three-body function, namely non-interacting electrons scattering 
against random impurities. The Hamiltonian representing impurity scattering is quadratic, 
and hence Wick's theorem is applicable, which means that the expectation value can be 
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factorized into pairwise contractions, i.e. expressed in terms of Green functions. Impurity 
averaging, which is now imphcit in the expectation value, reintroduces correlations between 
the particles, which implies that one must introduce vertex functions. However, we do not 
allow impurity correlations between the two subsystems, i.e. we assume that (AiA2)imp = 
(Ai)iinp(A2)imp-l^ The particular choice for the impurity self-energy used in the calculation of 
the impurity averaged Green function fixes the choice of the vertex function; in what follows 
we use the self-consistent Born approximation for the self-energy, and the corresponding 
vertex function consists either of the ladder diagrams (Section |III Q) , or of the maximally 
crossed diagrams (Section |IIID|) . The form for impurity- A giving the dominant contribution 
is shown in Fig. ^.Ill We consider only uniform systems, and set the external wave vector to 
zero, Q = 0. We also denote fermionic complex frequencies by ikm in contrast to bosonic 
frequencies iun, and the external frequency Q. Thus we have 



2 1 

A"(q, q, i^n + iuJn, ii^n) = X ^ X ^" ^' ^^«) 

+}C{k,-q,ikm,i^n,-{i^^n + i^n)) , (24) 



where 



^{k-i Qi ikjyii i^ni "i^n) — ^(^; ^km}'y{k, k] ik^, ikm + ^^n)^(^) ikm "I" i^n) 
xT{k, k + q; ikm + i^n, ikm + i^n + i^n)Q{k + q, ikm + i^n + i^^n) 
xT{k + q, k]ikm + ii.^n + i^nyikm) ■ (25) 



The factor 2 in Eq. (|2^) comes from the spin sum. In Eq. (|25D , k"-^ /m is the current (vector) 
vertex function and F is the charge (scalar) vertex function. In labeling the variables in the 
vertex functions, we use the convention that the incoming momentum (frequency) is first 
variable, and the outgoing momentum (frequency) is the second variable. The second term 
in the square brackets corresponds to reversing the order of the two U12 lines; see Fig. ^. 
Eqs.(p^)-(|25D need to be analytically continued to the real axis, after which they can be used 
as a starting point for evaluating the transconductance in the weak and strong scattering 
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limits, respectively. It should be noted that Eq. (^) does not include all possible diagrams. 
An example of a diagram not included is shown in Fig. ^(c). 

A. Analytic continuation 

The summation over the fermion frequencies follows the standard prescription:^ the 
discrete sum is replaced by a contour integration, [3~^ J2ikn fi'^^n) = — (27ri)^^ / dznp{z)f{z). 
To evaluate the contour integral, one must pay attention to the branch cuts of the integrand, 
and in the case of Eq.(|2^) these occur at Im[z] = 0, — f^m, — — (first term), and at 
Im[z] = 0,-^111,^^111 (second term). After performing the contour integration, we must 
further set iQm + i^m — > + + i5, iVLrn VL + i5, and iojm ~^ uj — i6 (see Eq. (|2^) ). After 
some tedious, but straightforward, algebra one finds 

2 rOO rjf: 

k 

x[K{k,q,e,VL,uj) + K{k,-q,e,VL,-{to + VL))^ , (26) 

where 

K{k, q, e, fi, u) = 

G"'(fc, e + n)T++{k, k + q,e + n,e + uj + n)G''{k + q,e + uj + n) 
X {r++(fc + q,k,e + uj + Q, e)G''{k, e)7++(fc, fc; e, e + Q) 
-r+_(fc + q,k,e + uj + n, e)G''{k, e)7_+(fc, fc; e, e + fi)} 
+ 

Cik + q,e + uj)T^4k + g, fc, e + cu, e - n)G''{k, e - n) 
X |7_+(fc, k;e — ^l, e)G''{k, e)r+_(fc, k + q,e,e + oo) 
-7_„(fc, k;e-n, e)G''{k, e)r__(fc, k + q,e,e + uj)} 
+ 

G°'{k, e - Q - uj)^_+{k, k;e - Q - uj,e - uo)G'^{k, e- u) 
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X |r++(fe, k + q,e — u;, e)G^(k + q, e)r+_(fe + q,k,e,e — Q — u) 

-r+_(fc, k + q,e-uj, e)G''{k + q, e)r__(fc + g, fe, e, e - 1^ - tu)} . (27) 

Here the subscripts ±± indicate the signs of the (infinitesimal) imaginary parts of the vertex 
functions' frequency arguments. This result is still quite general, and can be evaluated within 
different levels of approximation, of which we shall illustrate three special cases. 

B. The Boltzmann limit {ut > 1 and/or Dq^r > 1) 

In the weak scattering limit, we can neglect the charge vertex corrections, since F differs 
from unity only in a small region, where cut and Dq^r are small. Here r is the life-time due 
to impurity scattering and D is the diffusion constant. In the dc and weak-scattering (or 
"Boltzmann") limit, Ab becomes 

A§(q, q-u + i5,u-i5)^—Y. — n^{e)k- 

mu J-oo 2m 

x{KB{k,q,e,n = 0,lu) + KB{k, -q, e,n = 0,-lu)} , (28) 

where 

KB{k,q,e,Q = 0,u) ^ 

G'{k, e)G'{k + g, e + uj)[G'{k, e)7^+(fe, fe; e, e) - G\k, e)7^+(fe, fe; e, e) 
+G\k + g, e + uj)G\k, e)[7^?+(/e, k- e, t)G'\k, e) - 7^?_(fe, k- e, e)G\k, e)] 
+G'"(fe, e - a;)7?+(fe, fe; e - a;, e - uj)G''{k, e - uj)[G''{k + q,e)- G^ik + q, e)] 
^ 7^+(fc, k] e^, tj^){-G\k, e)G\k, e)[G'{k + q, e + cu) - G\k + g, e + uj)] 
+G''{k, € - uj)G\k, e - uj)[G'{k + q, e) - G\k + q, e)]} . (29) 

In writing the approximate equation, we used the fact that in the Boltzmann limit, (i) the 
terms of the form G^G^ or G°'G"' (with equal momentum and frequency arguments) are 
smaller by a factor of l/rEp <ti 1 than terms of the type G^'G"", and can hence be neglected 
in the present level of approximation, and (ii) G'^{k,uj)G"'{k,ijj) — TA{k,ijj), where A{k,ijj) 
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is the spectral function and l/2r = — ImS^(fc, ej^). A{k, u) is sharply peaked around e = e^, 
and hence the energy arguments in 7_+ can be replaced by ey^. 

In the Boltzmann limit of the current vertex function 7^ is well-knownS 

7^,(fc,fc;e;,,e;,) = ^. (30) 

Using G^{k, to) - G^ik, u) = -iA{k, u), G^{k, uj)G%k, u) = rA(fc, 00) and Eq. (H), Eq.® 
thus simplifies to 

A%{q, q;uj + iS,uj-iS) ^ V / — nF(e)fc°rtr(fc) 

k 

X [A{k, e) {A{k -q,e-uj) + A{k + q,e + uj)} 
-A{k, e + Lj)A{k - g, e) - A{k, e - uj)A{k + q, e)] 

X [{k + g)"rt,(fc + q) - + q, e + cu)A(fe, e) . (31) 

The Boltzmann limit is recovered by using free Green functions, which implies that the 
spectral functions reduce to 5-functions. We find 

Ag(q,q;c^ + z5,cu-z5) = — F"(qr,cu), (32) 

m 

where the transport polarization is given by 

F"{q,uj) = — Im^ ^ ^ [(fc + qrrUk + q) - fc"n,(fc)] . (33) 



Here rtr determines the in-plane conductivity, an = e^njrtr,i/mj. When (B^p is inserted into 
the expression for transconductivity, Eq. (P^j) , we obtain 



P2i = -^^ = -;^-E r 1^ ■ .^J Uu{q,u;)fF,{q,u;)F,{q,u;), (34) 

o"iiO"22 2nin2 ^ J-00 27r smh [pco'/2J 

Several comments are now in order. Without an applied magnetic field, the transconduc- 
tivity and consequently also the transresistivity are diagonal in the Cartesian coordinates, 
and we have suppressed the {a/5} indices in (0). For constant r's transport polarization is 
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related to the (bare) RPA polarization function, F°'{q,uj) = q'^lmxoiQ,^^)- In this limit Eq. 
(p^) reproduces the standard result for transresistivity, see, e.g. Refs. Since the above 
derivation ignores all higher order and/or quantum mechanical processes, it is not surprising 
that one can derive (^) directly from the Boltzmann equation.0 We also emphasize that in 
general the drag rate, or the transresistivity, cannot be expressed in terms of the polarization 
function; rather, one must use the more general object F"' defined above. 



C. Diffusive limit (ojt < 1, Dq'^T < 1) 

In this section we evaluate the transconductivity, in the weak scattering limit and in 
the diffusive limit {Dq^r < 1 and ujr < 1), including vertex corrections. Specifically, we 
consider momentum independent relaxation times, in which case t^^. = r, and include vertex 
corrections due to ladder diagrams. Then, we have 7 = 1 for all ± combinations, and the 
charge vertex is given byS 

m k + q,z„ z,) = 0[-lmiz,)lmiz,)] ^ ^^[lni(^,)Im(^2)] , (35) 

rjL'g^ - t{zi - 2;2)sgn[lm(2;i - Z2)\f 

where 6 is the step function. It is now straightforward to use ( P5| ) in the i^-function (P?!). 
Including only terms which involve G^C^ (with same arguments), and introducing a short- 
hand notation 

T^{q,uj) = [r{Dq^±iuj)]-\ (36) 

allows us to write 

Ki^{k, ±q, e,Q = 0, ±uj) = G'\k, e)G%k, e) 

X '-G''{k ±q,e± uj)T^{q, u) + G"(fc ±q,e± uj)T^{q, u) 

+G''{k, e T uj)G%k, e T uj) [G''{k ± q, e)r^(g, u) - G\k ± q, e)r±(g, u)] . (37) 

The triangle function in the ladder approximation, Ag, then becomes 

2t f°° de 

Al(9, q,uj + iS,uj -iS) = ^ fc" / — [nF(e + uj) - raF(e)] 

X {2Im[r-(g, cj)G'"(fc + g, e + uj)]A{k, e) - 2Im[r+(g, cj)G"(fc - q, e)]A{k, e + uj)} . (38) 
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We observe that the constant-r Boltzmann result is readily recovered from (^8]) by replacing 
r's by unity. A generalization to energy-dependent scattering rates is straightforward, but 
we do not reproduce the cumbersome results here. 

The next task is to establish a connection between the dressed polarization function 
x{q,uj) and the triangle function Al. In Appendix C we show that 

Imx(q, ^) = E / |^hF(e + uj)~ nF(e)]r-(g, u;)G^(/c + q,e + u)G\k, e)} . (39) 
k 

Comparison of (^) and (^) reveals some similarity, but clearly few more steps are required. 
We complete the connection by making a few observations. First, we express the spectral 
functions in (^) in terms of the retarded and advanced functions. The resulting integrals 
can be grouped in two classes: (i) integrals involving products of the type G""G", and (ii) 
integrals involving products of the form or G^G". We have earlier argued that type- 

(ii) integrals can be neglected in comparison with type-(i) integrals, if the momentum and 
frequency variables are equal. We now state that, in the present weak-scattering limit, this 
criterion applies also to functions with momentum variables and frequency variables which 
differ by less than Dt^^^"^ and (the diffusive limit). A proof for this statement is given in 
Appendix 0. Thus, keeping only the G^G"-terms in ( pHD allows us to express the quantity 
in curly brackets as 

{■ ■ ■} ^ - fr-(g, uj)G'\k + q, e + uo)G\k, e) + c.c 



+ 



r-(g, uj)G\k, e + uj)G''{k - q, e) + c.c 



In the above analysis the second term can be made to coincide with the first one by shifting 
the summation variable k ^ k -\- q; however when doing this the prefactor fc" in (pHf ) 
generates an extra q". This is exactly what is needed to give the required result, 

2rq" 

A^(qr, q,uj + i6,uj- i6) = -^lmx{q, uj) . (40) 

The above analysis shows the equivalence of the triangle function and the polarization 
function in the small q and u limit, confirming the result obtained by Zheng and MacDonaldi 
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with a different method. As observed by these authors, in the high-mobihty samples studied 
so far the replacement Xo X does not appear to be important; however, in dirtier samples 
the consequences of the vertex corrections {i.e. full x) i^^iy well become detectable. 



D. Weak localization correction to Coulomb drag 

In the previous sections we included the leading order impurity scattering diagrams which 
gave us the Boltzmann equation result for the case of weak scattering, and showed how the 
bare polarization function in a certain parameter range must be replaced by the dressed 
polarization function. Here we develop the analysis further and calculate the quantum cor- 
rection associated with weak localization (the basic physics of weak localization is reviewed, 
e.g. by Lee and Ramakrishnanil) . The corrections will be of the order of l/ikptj ^ 1, 
where ^ = vpT is the elastic mean free path. 

In Fig. ^ we display the different types of crossed diagrams that exist for the function A. 
The maximally crossed one is the one shown in Fig. ^(c). This diagram is, however, smaller 
than the one showed in Fig. ^(a), because of the restricted phase space. The two Green func- 
tions attached to the current vertex in the diagram in Fig. ^(a) have the same arguments 
because in the limit {Q = 0,Q = 0) the current vertex leaves the momenta and energies of 
the entering and leaving Green functions unchanged. Therefore there is the possibility of 
two overlapping spectral functions giving a overall factor of r. This does not happen for the 
diagram in Fig. §(c). Neither does the diagram in Fig. ^(b) lead to overlapping spectral func- 
tions except in all very small region of q, uj space, where q and uj are the incoming quantities 
at the charge vertices. Since we are integrating over g, uo the (logarithmic) singularity caused 
by the maximally crossed diagrams becomes regularized. In other words, the contribution 
from this diagram is small for the same reasons that the dressing of the charge vertexes, 
discussed in Sec. p.11 C| and can be neglected for experimentally relevant parameters.^. We 



therefore conclude that diagrams of the Fig. ^(a)-type dominate the quantum correction to 
the drag rate. 
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The leading quantum correction is given as the sum of the maximally crossed diagrams, 
the Cooperon. The resulting vertex function describing the weak localization correction, 
7^^", obeys 

WL 1 k ■ k' , , , 

uy [k'Y 

(41) 

where the Cooperon is given by0 

C{k,k]ikm,ikm + i^n) = - tt^ r, r, ; .„ . , (42) 

ivrpr 1 - ikm, ikm + 

where 

1 1 

C(fc, fc'; z/cm, ikm + «^^n) = 7^ ^(P + + iQn)Q{p, ikm) , (43) 

/vrpr z/ p 

where Q = k + k'. 

In order to evaluate the (^-function, we make use of the fact that the weak localization 
divergence occurs for small Q. With this in mind we replace e{k + Q) by e{k) +vfQ cos{6). 
Then we can integrate over k' in Eq. (^Sf ). For small DQ^ and Qr we obtain 

, if (fcm + ^^n)fcm < 

C(fc, fc ;«A;„,zA;^ + = < , (44) 

1^ 1 - |f]„|r - DQ'^T , if {km + ^n)km > 

where the diffusion constant is defined as D = VpT/2. This expression is only valid for 
Qi < 1, therefore the upper hmit of the integral in Eq. (|43|) has to be cut off by 

Next we perform the analytic continuations that are needed for the evaluation of the 
function A. Since the analytic continuation Un ^ ±i6 leads to \u!n\ l^iuj, we obtain for 
the ^-function 

C++ = C- = 0, 

C-+(fc, k'- e,e + n) = C_+(fc, k'; e-n,e) (45) 
= l + inT- DQ^T. 
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The Cooperon that enters the expression for Ky/i-, Eq.(^), then acquires the famihar form 

1 1 

" 2^^^nVDQ^' ^^^^ 
After integration over Q we find the weak locahzation vertex function, 

7^^(fc, k-e-VL,e)= -/"^^{k, k;e,e + Q) 

^ Aik, In(fir) = A(A;, e)r/WL(fi)/2r. (47) 

TlKFtT 

Here rf^^ is the ratio between the quantum correction and the classical conductivity: 
7]^^{Q) = S(t{uj)/(Jq^ The combinations 7^ and 7^ are both zero. We now get for 
the i^'-function 

Ky^Lik, q, e, Q, u) ^ -17]"^"^ {-[A{k , e)]^A{k + q,e + u;) + [A{k, e - uo)fA{k + q, e)} (48) 

Using that A^ tA/2 for large r we can express the weak localization correction 5 Awl in 
terms of the response functions as 

5^'^^{q,q;uj + i5,u-i5,^l) = ' ^ g"Imxo(g,^)=r/'^^(^)Ag(q,q;^ + 1(5,^-1,5), (49) 

which immediately leads to the conclusion that to leading order the weak localization cor- 
rection to transconductance is 

5aTm = [r/r(fi) + v^'-m^- (50) 
Consequently, the transresistance is unaffected by the weak localization correction because 



_ ^21(1 + ^1+^2) ^ . . 

f^ii(l + ^1)^^22(1 +^72) 



Weak localization is strongly affected by external magnetic fields. The formalism presented 
above can be extended to include magnetic fields; in particular, the topology of all diagrams 
remains unaltered. 
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IV. FINITE FREQUENCY RESPONSE AT T = 



A. General expression 



For finite frequencies and finite temperatures the analysis becomes considerably more 
complicated and for simplicity we therefore restrict ourselves to study the finite frequencies 
case at zero temperature. Furthermore, we will consider a model where all vertex functions 
can be replaced by unity, i.e. a system with short range impurity potentials and a system 
not in the diffusive limit. The transconductivity in terms of the time-ordered current-current 
correlation function is 



a"^(Q,n) = ^ (iRe [uf{Q,n)] - sgn(f])Im [uf{Q,n) 



(52) 



The time-ordered current-current correlation function is now written in terms of time-ordered 
Green functions as 



uf{Q = 0,Q) 
where 



[^-E\Ui2{q)fAt{q,q;n + u;,uj)A^{-q,-q;-n-uj,-uj), (53) 



Af (g, q-Q + uj,uj)= f^-Y, v^{k)G\k, + ln)G\k, - ^Q) 



X 



G\k + q, cui + cu) + G\k -q,uji- u] 



(54) 



It is straightforward to see that Af {q , q; u , u) = 0, which can be shown along the same 
lines as in the last part of Appendix 0. Furthermore we can show that A{q, q;uj,uj + Q) is 
proportional to Q, and consequently the transconductivity vanishes at zero frequency and 
zero temperature in agreement with Eq. (p2]). 



B. Clean system 



To evaluate the zero-temperature correlation function it is useful to decompose the 
time-ordered Green function into advanced and retarded parts according to G^{q,uj) = 
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Q{u!)G^{q,u!) + Q{—uj)G°'{q,uj). In the case of a clean system, the decomposition can also 
be carried out in momentum space as G'^{q, oj) = 0(|q| — kp)G'''{q, u) + Q{kp — \q\)G°'{q, u), 
which is actually more convenient since it leaves the frequency integrals unrestricted. Conse- 
quently, using the momentum space decomposition, we can carry out the conventional pole- 
position analysis and find that most of the terms arising from ( ^5] ) vanish. The remaining 
non-zero terms are most conveniently evaluated using the frequency-space decomposition, 
which yields for i7 > 

4-7 

nr(Q = o,fi) = -- y: \Ui2{q)\Wik)vM 

q,p,k 

^ J-n/2 2n J-n/2 27i Ju.^ 27i ^ ' ^ 2 ' ^ ' 2 ' 
xG' ik- q, ^1 - ^)G"-(p, u;2 + l^)G^ip, UJ2 - ^^)G''ip - ^2 - uj) 

+ / —G{k,ui + -Q)G{k,uji--Q) 

xG'ik- q, iu, - cu)G"-(p, UJ2 + ^n)G''ip, UJ2 - ^n)G''{p + q,iU2 + uj) 

+ / , IT —G'{k,uJi + -Q)G''{k,u,--Q) 

J-n/2 271 J-n/2 271 Jwi 27[ 2 2 

xG%k- q, uJi - uj)G'\p, UJ2 + ln)G''ip, t02 - ^l^)G"^(p + q,i02 + co)]. (55) 

Since all frequency integrations run at most over an interval of length Q, the end result 
is proportional to Q^. We can furthermore show that Il^^{Q,Q) = n^„(Q, — fi), so that 
Il'^^{Q,Q) ~ Using the fluctuation-dissipation theorem we find Re 11^ ~ \Q\^ and 

Im W ^ Q'^. Thus, to a leading order in r, the real part of the transconductivity is propor- 
tional to Q"^, and its imaginary part is proportional to sgn(f2)f2^. 

C. Disordered systems 

For disordered systems we can only use the frequency-space decomposition, and conse- 
quently the pole-position analysis is not quite as powerful as for clean systems. We can, 
however, determine the leading corrections by regarding perturbation and using 

a Taylor expansion of the type G^ {uJi)G'^ {uJi)G'^ {ui — lj) = GQ{LJi)GQ{iJi)Gl{LJi — u) + 
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{i l2T){d I dui)[GQ{ui)GQ{ui)GQ{uJi — oj)]. After the expansion, all propagators are given by 
clean system Green functions, and we can easily do the integral over one of the frequency 
arguments (in the example over uji). The resulting term is identical with a term that is 
encountered in the evaluation of an auxiliary time-ordered function for a clean system. The 
auxiliary function can be analyzed also by means of the momentum-space decomposition 
(since all propagators are given by Gq), and the frequency dependence of the various terms 
can be obtained in a manner similar to what we did in the clean system case. Carrying out 
the analysis for all terms arising from (|53|) , we find that the leading order corrections to 11'^ 
are of the form fi^/r, and we obtain 

n,(Q,fi) = Fo(Q)|^]|3 + Fl(Q)^]V^. (56) 

Thus, we finally have 

fi, |r]|<l/r 

a2i(Q = 0,r])~ (57) 
[ n\ 1/r < \n\ ^tp/h 

The constant -Fo(O) can be evaluated approximately be keeping only the most important 
terms. Taking only terms that are leading order in r into account, we find 




V. DISCUSSION AND CONCLUSIONS 

In this paper, we have presented a fully microscopic theory of the Coulomb drag, based on 
the Kubo formalism. We have used the finite-temperature formalism to obtain expressions 
for the dc drag, and the zero-temperature formalism to obtain finite frequency results. We 
have chosen to present only formal results here, deferring the presentation of experimental 
consequences of these results to another publication.0 

We calculate the transconductivity ai2 using an order by order expansion in the inter- 
layer interaction U{q). Assuming no correlations between the impurities between the two 
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layers, we find an exact relation between the first order result a^2{Qj^)j the subsystem 
conductivities. The result also indicates that in a uniform system, the dc transconductivity 
vanishes to first order. 

(2) 

To second order, we write a formal result for the transconductivity crj2 terms of 
the A(q, q'; cu, to'')-functions, which are the thermal-averaged (jpp) correlation functions of 
the individual subsystems. In evaluating A under various circumstances, we find (i) for 
constant intralayer elastic scattering rates, we duplicate in the limit 1/t —>■ results obtained 
earlier using the Boltzmann equation, and the memory functional method in the diffusive 
limit; (ii) for energy-dependent intralayer elastic scattering rates, however, the g"Im[x] 
must be replaced by another quantity F°'{q,uj), which we call the transport polarizability. 
The energy-dependent result is due (from the Boltzmann equation point of view) to the 
the fact that the perturbed distribution function on application of the electric field for 
energy-dependent elastic scattering rates is not a drifted Fermi-Dirac.0 However, intralayer 
electron-electron interactions tend to relax the distribution function back to a drifted Fermi- 
Dirac, and hence the larger the intralayer e-e interactions are, the closer the F°'{q,uj) will 
be to g°Im[xo] in Eq. (0).^ 

We have calculated the weak-localization correction to the second order transconductiv- 
ity, and find that ^cr^-^/cr^i = 5crn^/o"ii + ^cr^^/cr22, which implies that, to lowest order in 
(kpi)'^, the transresistivity p2i is unaffected by the weak-localization corrections. 

The zero-temperature formalism indicates that the dc-drag vanishes in an open system. 
This result is reproduced in the T ^ of the finite temperature formalism, since the di^n^{uj) 
term in Eq. (|22D vanishes as T ^ 0, and one can show that A is linear in u; as — > 0. This 
statement is valid for open systems in which are connected to dissipative leads, and not for 
closed systems.lli 

For finite frequencies we have evaluated the leading contribution to the transconductivity 
at zero temperature and found that in the clean limit o"2i ~ fi^. Including disorder we showed 
that frequency in this expression is replaced by 1/r and a2i ^ Vt/r. 
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The formalism that we have developed in this paper can be applied to many different 
physical realizations of coupled electron systems. It is thus straightforward to extend the 
calculation to include magnetic field and work is in progress in this direction.B The present 
formalism also forms a useful starting point for the study of e.g. higher order intralayer inter- 
actions, phonon mediated intralayer interaction, and correlations caused by strong interlayer 
electron-electron interactions. 
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APPENDIX A: BRANCH CUTS FOR THE THREE BODY CORRELATION 

FUNCTION 

Consider the Fourier transform in time according to Eq. (pIsD (we set r" = since A is 
only a function of r — t" and r' — t") 

A{inn + iujn, i(^n) = 1^ dr r dr' e'^"^-^''^"^' A{xx' x" ; tt') 
Jo Jo 

= dTe'"""^ dr' e''^"^' {J{xT)p{x'r')p{x")) + J^^ dr' e'^"^' {p{x't')J{xt)p{x"))^ . (Al) 

We now insert the identity 1 = J2n where {|^)} is a set of eigenstates for the Hamilto- 

nian, between the operators in Eq. (|AT|). After performing the two imaginary time integrals 
and some algebra we obtain 

1 



A(in„ + = e''^^ {k\J{x)\m){m\p{x')\l){l\p{x")\k)- 



X 



kml ^ + Em-Ei 

.iflri + i^^n + Ek-Ei^ ^ lVL^ + Ek- Err,,^ H 
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+ {m\J{x)\k){l\p{x')\m){k\p{x")\l) 



iun + El- E, 



m 



X 



[iQn + iujn + El - Ek ' ^ tQ^ + E^-Ek^'' " ' e ^5 ")}) (^2) 



From this expression we can read of the branch cuts of A{iQn + i^n-,i^n)i which is used 
when doing the Matsubara sum over iojn in Eq. ( pOa| ). 

APPENDIX B: PROOF OF A(++) = 
We consider the quantity 

A{x',inn + iuJn,iuJn)= f dr' rrfre*^"V^""'(T(J(r)p(3;V)p(0))), (Bl) 

Jo Jo 



where J = J dxj^(x). The Fourier of A{x', ifln + i^^n, i^n) with respect to x' is A{q, q; ifln + 

We need the two combinations A(h — ) and A(++). It is clear that if we set = in 
expression above we lose the information necessary to evaluate A(H — ) and we can only get 
A(++) by the substitution iun u + i6. The dc-limit of A(++) can be safely evaluated 
by setting iQn = before the analytic continuation. Doing that we obtain 

A{x\tuJn,tuJn) = f3 dr' e''^"^' {T{J{r)p{x'r')p{0))). (B2) 

JO 

Now consider the density- density correlation function for the case when the Hamiltonian 
has been enlarged by a vector potential term 

H^H + A^J^. (B3) 

(Omitting a diamagnetic term which is not important for present argument.) If we view 
as a perturbation we can write the charge-charge correlation function as 



Xa{,x\ t') = (r(exp ^- ^ drA^J^ 



(r)Jp(xV)p(0))). (B4) 
From this expression it is seen that the function A can be obtained as 
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(B5) 



Since a constant vector potential can always be removed by a gauge transformation (this is 
however only strictly true for a system with open boundary conditions), xa can not depend 
on A, and hence we arrive at the conclusion that A = 0. 

As a specific example we now take the impurity averaged A(++)-function in the sim- 
plifying case where we can neglect the charge vertex corrections (i.e., not in the diffusive 
limit). Setting iQn = in Eq. (|25|) , we then have 

— /C(fc, q, ikm, iuJn) = [Gik, ifc^)]^ — l{k; ikm)Q{k + q, ik^ + iuJn) 
m m 

= ^TT- y{k + q, ikm + lUJn), 

ok"' 

where we have used the Ward identity: dj^G^^{k,ikm) = kj{k,ikm)/m. Integrating by 
parts, we obtain 

A"(qr, q, ; iu^, ^u;„) = t— — -^ Y1 / 7TrT?(^(^' ikm)G{k + g, ikm + iu;„) 



-Q{k, ikm)G{k - q, ik^ - iuJn)j = 0, (B6) 
which can be seen by shifting ikn and k in the term. 

APPENDIX C: IMAGINARY PART OF THE POLARIZABILITY IN THE 

DIFFUSIVE LIMIT 

In this appendix, we derive expressions for Im[x] when qvpT < 1 and ujt < 1; i.e, in the 
diffusive limit. The polarizability including the vertex correction F is given by 

2 

X{q, iuJn) = — ^ XI XI "^(^ + 9, i^n + iuJn)G{k, i€n)T{q; iSn + iuJn, i£n)- (CI) 

For qvpT < 1 and ujt < 1, F is given by Eq. (|35|) . Inserting this form of the vertex correction 
in Eq. ( |d|) , writing the sum as a contour integral and deforming the contour in the standard 
mannei0 yields, 
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+ £) - r-(g, cu)G"^(fc, e)]G'{k + q,e + oj)]. (C2) 

In the above equation, the analytic continuation iuj^ uj + iO~^ has been taken, and r^(g, u) 
is defined in Eq. (pG]). 

We write x = Xa + Xb, where Xa (Xfe) excludes (includes) the vertex corrections. The 
term Xa to lowest order in q and u is 
2 /■ r/r 

Xa(g, uj) = --Y.J ^n^ie) [G'-{k, e)G'{k + q,e + uj)- G\k + q, e - u)] 



J ^n',{e)lj:[G'ik,e)-G'^ik,e)] = -^ (C3) 



Note that Xa to lowest order in q and uj is purely real. 
The second part is 

2^ r de 



2 r (IF 

Xb{q, uj) = -r-(g, 00)- J2 TT-M^) [G''ik + q, e)G\k, e - u) - G%k, e)G\k + q,e + uj)]. 
u J 27ri 

(C4) 



To lowest order in q and to, this i: 

Xb{q,uj) = —-— — . C5 

a/i Dq'' — lu) 

Thus Xb has an imaginary component which, as shown in section pTl] , is related to the A 
when qvpT < 1 and ujt < 1. 

APPENDIX D: JUSTIFICATION FOR NEGLECTING G^G^- AND G^G^-TERMS 

In this appendix, we show that the terms in A^, involving products G^G^ and G'^G"' 
(denoted below as A^) can be neglected when compared to G'^G^ when EpT ^ 1, in the 
diffusive limit. 

Written in full, A^ is 

27 



2t 1 r dp 

A^(g, q; 00' + z6,u;' - z6) = ^ / 

m V , J Ztt 



2t1 ^ f de 

k 

r-{q, 00) \-G''{k, 6)G'\k + q,e + u;)- G^{k, -q,e-Lo) 



+G''{k,e-uj)G^{k + q,e) + G\k,e + uj)G\k-q,e) + c.c (Dl) 



Note that at = 0, this term is identically zero. Expanding in powers of uj gives 



A^(g, g; a;' + i5, J - ib) = Re 



mu , J In 



d[G^{k, e)G%k -q,e)- G^{k, e)G'-{k + q, e)] 



Re 



de 

-c^r- V / —n'p{e)k 



[^"(fe, e)G''{k -q,e)- G'{k, e)G'{k + q, e)] 



+ Oiu;^). (D2) 



Expanding Eq. (P2|) in powers of q, and assuming that the self-energy is small (since r is 
large) so that dG/dsj^ ^ G^, yields 



A^(g, q; u' + i5, oj' - i5) = Re 



u;q-T-{q,co)--Y. [ 



de , {k 



a\2 
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np[e) 



m 



[G"^(fc,£)3 + G"-(fc,£)^ 



+ 0{q\u') 



(D3) 



Integrals of e)^ over k do not diverge when t ^ 00 because the poles of the function 

are on the same half-plane (unlike integrals over G"'{k, e)G^{k, e), which go as r). Since n'p{e) 
is peaked around /x, one can estimate the magnitude of A" by replacing —n'p{e) ~ 6{e — 
which gives 



A'^{q,q';uj' + i6,uj' -i5) ^ Re 



r {q,uj[ 



erq u! 
2tt^Ep 



(D4) 



A comparison of this with A^ given in Eq. (|40|) shows that A" is smaller by a factor of 

l/{EpT). 
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impurity processes can safely be neglected. 

One should note that Fig. ^ does not include all possible diagrams. Examples of such 
higher-order diagrams are analyzed in Section [illlJ| . 

See, for example, Ref. 0, Section 7. I.e. 
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FIGURES 

FIG. 1. Diagrams corresponding to the current-current correlation function, to (a) first and 
(b) second order in the interlayer Coulomb interaction. The shaded triangles correspond to the As 
given in Eqs. (|^) and (|l8|), the dashed lines to the interaction, the dotted lines to the external 
current operators, and the arrowheads to the direction of momentum and energy transfer. 

FIG. 2. Diagrams which lead to the screened interlayer interaction within the random phase 
approximation. The "bubbles" are the bare polarizabilities of the subsystems, the thin wavy lines 
are the bare interactions, the thick lines are the screened interactions, and numbers indicate the 
subsystem. 

FIG. 3. The function A for the case in which vertex corrections are included at each of the 
individual charge and current vertices. Figure (a) shows the decomposition of A into diagrams 
with clockwise- and anticlockwise-moving Green functions, with the grey shaded areas indicating 
vertex corrections. Figure (b) shows one of these diagrams in greater detail. Here A;"7 is the current 
vertex, the F is the charge vertex, the dashed line is incoming momentum and frequency, the dotted 
lines are the interaction Ui2{q,uj), and the solid lines with arrows are the Green function. Normal 
momentum and energy conservation rules apply at the vertices. 

FIG. 4. The different types of crossed diagrams for the triangle diagram relevant for the weak 
localization correction to the transconductivity. Diagrams of the type in (a) give the leading order 
contribution. Dressing the charge vertexes as in (b) gives a smaller contribution for moderately 
clean samples for the same reason that allows us to neglect the vertex corrections of the charge 



vertexes, discussed in Sec. QIC. The diagram in (c), which cannot be included using vertex 



functions alone, has an even smaller phase space and can hence also be neglected. 
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